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Abstract 

We present a general procedure for calculating the exact partition function of an Ising 
model on a periodic chain in presence of magnetic field considering both open and closed 
boundary conditions . Using same procedure on a quasiperiodic (Fibonacci) chain we 
have established a recurrence relation among partition functions of different Fibonacci 
generations from n — th to (n + 6) — th . In the large N limit we find (2r + l)F n+ % = -F n -2 
; where r is the golden mean and F n stands for free energy/spin for the n — th generation 
. We have also studied chemical potential in both cases. 
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I Introduction 

In this paper we have calculated the exact partition function of an Ising model on a 
periodic chain in presence of magnetic field with open boundary conditions by splitting the 
transfer matrix into two particular non-commuting matrices. The result is quite nontrivial 
in contrast to the expression for the partition function of the closed one[l,2]. Since there 
has been an enormous amount of work on Ising model it is difficult to do justice to all 
references, however rigorous results and developments can be found in [1,3] and in the 
review article [4] We have studied chemical potentials for periodic chain in presence of 
magnetic field . Studying Ising model on a Fibonacci chain in presence of magnetic field 
is interesting since the discovery of quasicrystals [5]. Scaling forms of thermodynamic 
functions for such system have been studied using renormalisation group technique [6] 
through one step decimation . The ground state and thermodynamic properties of such a 
system have been studied using renormalization group technique [6,7] . Our formulation 
helps us to express the partition function of Ising model on Fibonacci chain in presence of 
magnetic field as a sum of partition functions of usual Ising open chains with coefficients 
containing Fibonacci symmetry. We have also studied some symmetry properties of the 
Fibonacci chain. Using a special symmetry property ("Mirror Symmetry") and the usual 
trace map relation [8,9] ( trace maps and invariants relating to two- letter substitution 
lattices have been studied in [10] and references therein ) we have established a recurrence 
relation among the partition functions of different Fibonacci generations. This includes all 
the partition functions starting from n — th up to (n + 6) — th generations. We observe that 
mirror symmetry is a characteristic property of each Fibonacci generation with n — th and 
(n + 6) — th generations having same topology [11]. Assuming that the free energy/spin 
for both open and closed chains for a particular generation are equal in the large N limit 



we have obtained a scaling property among the free energy/spin of two consecutive even 
genarations, using the above recurrence relation . In this particular case we have studied 
chemical potential in absence of magnetic field. 



II Exact partition function for open Ising chain with magnetic field 

The one dimensional Ising model consists of a chain of N spins Si — ±1 

;i — 1, 2, , N with nearest neighbour interactions e^+i. The 

Hamiltonian is given by: 

N-l N 

7~L = — X! e i,i+iSiSi+i — H^Si (1) 

i=l i=l 

For a uniform lattice e iji+ i = e,the partition function is given by: 

Z° N {T, H) = £ f( S ^ W(S 2 , S 3 ) f(S N _ u S N )f (S N , SJ (2) 

Si,S2,----,Sff=— 1 

with f(Si, S i+ i) = exp[(3eSiS i+1 + \(3H(Si + S i+1 )]; f (S N , Si) = [f(S N , 5i)] e=0 . Here 
the superscript o stands for the chain with open boundary condition. Therefore the parti- 
tion function (2) can be written in terms of transfer matrix as: 

Z° N (T,H) =TrP N - 1 P (3) 

where 

P = v^(l + £)ffi = vWl + ^) ( 4 ) 

P = [P}e=0 = (1 + A) (7i = ^(1 + X T ) (5) 

with r = exp(-2(3e) , A = (y e -pH 6 {]\ XT= ^^ H 6 Q ) and(T i = (i J 

The general formula of the partition function for even and odd number 
of spins (i.e., odd and even number of bonds) can be derived by using equations 
(3) and (4) as : 



and 



where 



Z° 2N (T, H) = r N ~Tr{l + Xl ){\ + x 2 )....(l + x 2N . 1 )(l + \ T ) (6) 
(T, H) = r N Tr(l + Xl )(l + x 2 ) (1 + x 2N )(l + X)a L (7) 



x 2 i+i = -, x 2i = — ; i = integer (8) 
r r 

The above equations show that A, A T are the signatures for the transfer matrices 
corresponding to bonds in odd and even positions. In the case of a chain with closed 
boundary condition the last factor in eqn.(2) is f(S^, Si) and consequently the partition 
function takes the form: 

Z C 2N (T, H) = r N Tr{\ + Xl )(l + x 2 ) (1 + x 2N ) (9) 



2 



Z 2N+1 = r N+1 *Tr(l + Xl ){\ + x 2 ) (1 + x 2N+1 )a 1 (10) 

Here the superscript c indicates closed chain. One can show by elementary calculation 
that eqns.(9) and (10) reduce to the well known form[12] 

Z C N {T,H) = \ N + +\ N _ (11) 

where 

A ± = r - 3 [cosh (f3H) ± ^(sinh 2 ^) + r 2 )] (12) 

are the eigenvalues of the transfer matrix P. The expression for the partition function 
in the case of an open chain with even number of spins can be derived from eqn.(6) as 
follows: 



Z° N (T, H) = r N - l iTr(l + Xl )(l + x 2 ) (1 + x 2N _ 1 )(l + A T ) 

= V^Z C 2N (T, H) + r N ^(l - t)Tt{1 + Xl )(l + x 2 ) (1 + X2N -i) 

= V^Z C 2N (T, H) + v^(l - r)Z 2 c (iV _i)(T, H) + r N ^{l - r) 

xTr(l + Xl )(l + x 2 ) (1 + x 2N _ 2 )x 2N ^ (13) 

The last term in the above expression can be written in terms of the eigenvalues of 
the transfer matrix P viz.A±. By following the method of induction : 



r N ~*(l - r)Tr{\ + ^(l + x 2 ) (1 + x 2JV _ 2 )x 2JV _ 1 

= (l-r)r^Acosh 2 (^)E(^-) (~ r 



} _ X J(N-1)_ X _2(N-1) 



4(1 _ r y- 1 cosh" \/3 'H)— — —3 (14) 

A_l — A_ 



So eqn.(13) becomes : 



Z° N {T,H) = V^ZU^^ + V^il-^Z^^H) 

\2{N-1) _ X J(N-1) 

+4(1 - r)r~ cosh 2 {(3H) x — — — 2 — 2 (15) 

A_|_ — A_ 

Similarly the expression (7) for the open chain partition function with odd number of 
spins takes the form: 



Z 2 n+i 

(T, H) = V^Z C 2N+1 + 4(1 - r)cosh ((3H) 

\ 2N \ 2N 

x A t (is) 

A + - A + 

It can be shown easily that the free energy per spin calculated from equations (11) 
, (15) and (16) are all same. Thus the forms of the thermodynamic functions are same 
for both open and closed chains . However as the expressions of the partition function 
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are different for closed and open chain it will be highly instructive to study the chemical 
potential of such systems. 

Ila Chemical Potential 

The chemical potential for an open chain with N number of spins is given by: 

_ d logZ° N _ Z° N+1 
fiN ~ dN ~ L ° 9 Z° N U7j 

The expression for the chemical potential for a system with even number of spin ' say 2iV 
; is obtained by adding a spin to the spin to the system and substituting from equations 
(15) and (16) in the above relation (17) in the large N limit: 



Z° N+l _ 2 v^c) + 4(1 - r)cosh(/3H) x 



"even ~ 70 ~ \ u c) * 



Z° 2N ^F(v c ) 2 + ^F(l-r)+4r-(l-r)cosh 2 f3H 



z c 

where v c = e Mc = y C +1 = A+ , u c being the chemical potential for the closed chain. For 

N 



H = 0, A + = 2cosh f3e and consequently expression (18) becomes 

v° even = 2cosh Pe (19) 
Similarly for a system with odd number of spins we can write 

= < dd = ^ = x (20) 

Zj 2N-\ u even 

Therefore, v° even x v° odd = (u c ) 2 which implies 

f^even + l^odd = ^^c (21) 

From the above expressions we conclude that the chemical potentials ^c-,l^ even ^° dd are a ^ 
different. However, for H = 

A*c = P°even = Vodd = l °9 ( 2 cosh (3e) (22) 

which shows that the chemical energies for three different conditions of the chain men- 
tioned above are degenerate and it is removed when magnetic field is applied [equation 
(21)]. 

Ill Ising model on Fibonacci chain with magnetic field 

A Fibonacci chain can be inflated by two bonds L{large) and S(small) by the inflation 
rule L — > LS, S — ► L.The chain can be represented by the sequence: 

L — > LS — ► LSL — ► LSLLS — »■ LSLLSLSL — ► (23) 

In this case the interaction strengths in the Hamiltonian (1) e iji+1 = e for long bonds 
and e^j+i = e for the short ones where the bonds are arranged according to the Fibonacci 
sequence (23). The corresponding partition function of the nth generation Fibonacci chain 
having iV spins with N — 1 bonds is given by: 

Z° N (F) = TrPPPPP P (24) 
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where for long bonds the transfer matrix P is given by eqn.(4) and for short bonds the 
transfer matrix P is given by eqn.(4) with r replaced by f = r| e= -. Henceforth Z^(F) and 
Zpf(I) will represent partition functions for Ising models on an open Fibonacci chain and 
on an open regular lattice respectively. The expressions for the partition functions with 
odd and even number of bonds take the same forms as shown in eqns.(6) and (7) with x^s 
given in eqn.(8) for long bonds whereas for short bonds we replace r by f in eqn.(8).The 
explicit expressions for the partition functions for open and closed chains are: 

JVr Nq _ 

Z° 2N {F) = r—f-Tr{\ + ^(l + x 2 ) (1 + x 2N ^)(l + A ) (25) 

JVr Nq 

Z°2N+i{F) = r—f-Tr{l + ^(l + x 2 ) (1 + x 2N ){\ + X)a, (26) 

Note that the subscript 2N in the left hand side of eq.(25) stands for number of spins 
so that number of bonds is 2N — 1 which is odd. In eq.(26) number of spins is 2N + 1 and 
number of bonds is 2N which is even. 

Similarly for closed chain eqs. (24), (25) and (26) take the following forms: 

Z C N (F) = TrPPPPP P (27) 

JVr +1 No 

Z C 2N (F) = r—f-Tr{\ + x x ){\ + x 2 ) (1 + x 2iV _ 1 )(l + x 2N ) (28) 

JVr+l JVg 

Z 2 N+l(F) — r 2 f 2 Tr{l + X\){1 + X 2 ) (1 + X2Ar)(l + x 2 n+i)cti (29) 

where Nl,Ns are number of long and short bonds in a particular sequence. 
Now P is related to P through the following equation: 

P = y/f(l - 3)<Ti + ^P (30) 

Using eqn.(30) in eqns.(25) and (26) the Fibonacci partition function for any genera- 
tion can be written in terms of open Ising partition functions as follows: 

N 

Z° 2N (F) = h (e, e) + J2h2i(e, e)Z° 2i (I) (31) 
i=i 

N 

Z° 2 n-i(F) = lo(e, e) + ^w-i^ t)Z° 2i _ x {I) (32) 

i=i 

where Z 2i (I) and Z 2i _ x {F) are given by eqns.(15) and (16) respectively. We observe that 
the quasiperiodic nature of the Fibonacci chain is encoded in the functions h(e, e) and 
Z(e, e). Though for small generations these functions can be derived exactly still we could 
not find out their general forms. 

Ilia Recurrence relation among partition functions 

To circumvent the above difficulty we study the recurrence relations among the par- 
tition functions of different Fibonacci generations. A servey of different Fibonacci gener- 
ations depicted by eqn.(23) shows a symmetric pattern in terms of the number of bonds 
,viz., 
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Pi : P{odd) 
P 2 : PP{even) 
P 3 : PPP(odd) 
P 4 : PPPPP(odd) 
P 5 : PPPPPPPP(even) 
P 6 : PPPPPPPPPPPPP(odd) (33) 



and so on. 



The trace map relation [8,9] was introduced to study the spectrum of ID Schrodinger 
equation in a discontinious quasiperiodic potential. The use of tree map in different kind 
of substitution lattices have been studied in the review article [4] with many relevant 
references therein. 

Let P n _2 be the (n — 2)th Fibonacci generation with even number of bonds. This 
automatically ensures that the previous as well as the next two consecutive generations 
will have odd number of bonds . The recurrence relation for Fibonacci generations is 
given by: 



n- 



-lPn-2 (34) 

Now adding a term P„_ 2 Pn-3 in the above equation gives 

Pn + Pn-2Pn-3 = Pn-lPn-2 + ^-2^-3 (35) 

where P„_ 2 = Pet(P„_ 2 ).The following operations are applied sequentially on eqn.(35): 

substitute P n -i X Pi-i i n place of P„_3 on the right hand side and finally use Cayley- 
Hamilton theorem to get the usual trace map relation [8] on the Fibonacci lattice: 

TrP n = TrPn^TrPn-2 - D n _ 2 TrP n - 3 (36) 

The above equation will be necessary for calculating recurrence relation among dif- 
ferent Fibonacci generations. For this purpose we must understand symmetry properties 
of Fibonacci chain. Inspecting different generations of the Fibonacci chain it reveals that 
if the total number of bonds N of a particular generation is odd then there is a mirror 

reflection symmetry arround the (^j^Jth bond ;except the last two bonds. If the special 

bond arround which mirror symmetry occurs is a short (long) one the Fibonacci generation 
will have equal number of odd and even short (long) bonds . However if the total number 
of bonds N is even,the mirror reflection symmetry is arround a cluster of two successive 

long bonds at the ^y^t/i and (^j^jth positions of the chain . So "Mirror reflection 

symmetry" is a characteristic property of a Fibonacci chain. 

The nth and (n ± 3)th generations have the mirror reflection symmetry property 
arround the same kind of bond with last two bonds interchanged , while the (n ± 6)th 
generations are topologically same as the nth one. 

Using recurrence relation (34) we can write 

(37) 

Using Cayley-Hamilton theorem on the right hand side of eq.(37) we get: 
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Ai-2-P«-3 — {TrP n _ 2 )Pn-l ~ P n -lPn-\ (38) 

Multiplying eq.(38) by P from the right and taking trace we obtain: 

D n _ 2 Z c n _ 3 = (TrP n ^)Z c n ^ - Tr(P n ^P n ^P) (39) 
In a similar fashion we obtain : 

D n _ 2 Z° n _ 3 = {TrP n _ 2 )Z° n _ x - rr(P n _ 2 P„_ 1 P ) (40) 

The avove relations cleary show that one cannot obtain a recurrence relation among 
the partition functions through the trace map relation alone. To circumvent this difficulty 
we take recourse to mirror symmetry property. The expression P n - 2 Pn-\ i n eqns.(39) and 
(40) is similar to P n = P n -\P n -2 with last two bonds interchanged, i.e., both of them have 
the same mirror symmetric part fl n . Therefore eqns.(39) and (40) can be written as: 

D n _ 2 Z c n _ 3 = Z c n _ x {TrP n _ 2 ) - Tr(Q n PPP) (41) 



D n _ 2 Z° n _ 3 = Z° n _ x {TrP n - 2 ) - Tr(Q n PPP ) (42) 

The transfer matrix has the property that P = P T and P = P T .If such transfer 
matrices are arranged in a mirror symmetric fashion then the resulting matrix (Q n ) will 
have the following properties: 

i) Off diagonal elements are same i.e.,(u; n )i2 = (u) n )2i 

ii) Diagonal elements are not same but satisfy the condition: 
(vn)u(p,q) = (uj n ) 22 (q,p); where p = e pH ,q = e~ m . 

Thus the matrix Q n in eqns. (41) and (42) is of the form: 

q = f (<^n)ll (Vn)l2\ 

V(o; n )2i {u n ) 22 J 
Eqns. (41) and (42) can be written explictly in the following way: 



D n _ 2 Z c n _ 3 = Z c n _ 1 (TrP n - 2 )-rVf[Mii(y + — ) + M 12 (u + v + ^^-) 

r r 

+M 22 (x+^)} (43) 



D n _ 2 Z° n _ 3 = Z° n _ x {TrP n - 2 ) 
+(w n )22(x + qv)] 

where we have used 

p P = ^f( u y ) 

\X V J 

with 

x = 2 + i y = 2 + l u = i + 4 and v = 1 + 4. 
Eliminating TrP n _ 2 from eqns. (43) and (44) we have: 



- Vrr[(uj n ) u(y + pu) + {uj n ) 12 (u + v + px + qy) 



(44) 
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yVn^+puV'n-! xVn^ + qvV'n-! 

^ n + a,. , „,M/ . _l a™. _l -^22 + ^ 12 



(it + v)K-i + (pz + ?y) V"„-i (« + ^)K-i + ipx + qy)V' n -i 

D 7° 7 C 7 C 7° 

rv^ (u + v)V n _ 1 + (px + qy)V' n _ 1 { 

where 

K-i = ^ " ^-i and ^ n -i = ^ - 

To solve for (u; n )ii,(o> n )i2 and (u; n ) 2 2 we need another two equations. These equations 
are obtained from the usual formulae: 

Z c n = Tr(Q n PPP), Z° n = Tr(Q n PPP ). 

The explicit forms of these two relations are: 

x + ^ y + — 

'^n)ll H ; ; r ax+vv (u n )22 + (^n)l2 



K 1 (46) 



x + pu y + qv 

Wn)n H ; ; ; KJ22 + (^n)i2 



w + v + qx + py u + v + qx + py 

Z° 1 



(47) 



y/rr u + v + qx + py 
By elementary calculation one obtains: 

V — V Jl 

K)n(z,y) = p — 77 — 7 — x — 7^ x A n _i 

' ' x^ZX-^Jif) (48) 



y/rr r n _iA'„_i - r'„_iA n _i n a/t 



/ s / \ A n _i — A' n _i -Dn-2 a 

K) 22 (x,y) = x — -=■ x A n _x 

1 n _iA'„_i - 1 ' n _iA„_i rVr 

+ 4=T^V ^ A X i A ^ Z n K ' - A 'n~i%K) (49) 

Vrr r n _!A' n _i - r' n _iA n _! 



and 



(co n ) 12 (x,y) 



(qx + py)(l - }) 



xV n -(l- -)(pu(u n ) n + qv(u n ) 22 ) (50) 



'rr 



where 



An-i(x,y) 



yV n -i +puV' n ^i 



X + 



pu 



(u + v)V n ^ + (px + qy)V' n ^ u + v + ^ 



(51) 
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A' n _i(a;,y) = 
r n _i(x,y) = 

r / n-i(aj,y) = 



x + pu 



(u + v )K_i + (px + u + v + qx + py 



(« + + (px + qy)V' n ^ u + v + ^ 

xV n -i + qvV'n-! y + qv 



(u + v)V n -i + (px + gy)y n _i u + v + qx + py 

ryo rye <J C <Jo 

Zj n-\ Zj n-3 ~ ^n-l^nS 



K'(x,y) = 



(it + u) V n _i + {px + gy)y„_i 
1 



U + V + 



qx+py 



u + v + qx + py 
Eliminating D n _ 2 from eqns.(43) and (44) we obtain: 



(52) 
(53) 
(54) 
(55) 
(56) 
(57) 



TrR 



n-2 



70 yc 7c 70 

+ (o^ + P' 

1 



, v. 


-3 


xy y 
v n- 


-3 


V' n 


-3 


xy y 

' n- 


-3 



(u; n ) 22 (x,y) 



rr(l - \)(qx + py) 



u + v) + (px + qy) V r n 3 ) 1 , 



v; 



n— 3 



where 



u+v 



(x-y)(q-p) 
py+qx 

_ (x-yUg-p) 



x pu 

x 



(58) 



py+qx 11 r- xy py+qx 

In general nth and (n±2)th generations have same arrangement of the last two bonds 
appart from their respective mirror symmetric parts. That is why TrP n and TrP n ± 2 will 
have similar expressions. Since we have assumed P n = Q n PP it follows from the expression 
(33) that P n _3 = Q n ^ 3 PP. Therefore proceeding in a similar way as above we get: 



rvr x V r 



n— 4 



Zo yc 



n-2 ZJ n-4 



J rt-2 ZJ n— 4 



V'n-4 



yX y 



n— 4 



n— 4 



n-4 



rr(l - -){px + qy) 



u + v) + (qx + py) 



(w„_i)n(y,x) 

(Wn-l)ll(y,^) 



n— 4 



(59) 



Using eqns. (58) and (59) and similar expressions for TrP n , TrP n _\ in the trace 
map relation (36) we obtain the following recurrence relation among partition functions 
of different Fibonacci generations as: 
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70 7c 



7c 70 



a 



X- 



Ki-3 



rr(l - 



+ Pxy-^-^J (uj n +2)u(x, y) + (a' xy + Z?'^- 
u + v + (px + <?y) ^-^) V n+2 

Vn-l / 



n-1 



n-l 



u; n + 2 ) 22 (a;,y) 
x V n -2 



Z°Z, 



n-2 



ztzz 



n-2 



70 7c 



n— 3 



7c 70 
1 



a 



yy.i; 



+ Pyx~ry~ (uj n+l ) u (y, x) + (a' yx + (3' — —) (a; n+ i) 22 (y, a) 

*n-2 / \ Vn-2 / 



rr(l - i)(px + qry) 

X2/ "I" P xy T.r 

Vn-3 

A1-2K1-4 



70 7c 7c 70 

Z 'n-2 Z 'n-4 — Z 'n-2 Z 'n-4 



a 



It + f + + p|/) 

(u; n ) 22 (a:,y) - 



n-2 



n-2 



v r 



n+1 



X 



y.r 



)/:/■ 



n— 4 



(w„_i)ii(j/,x) + (a' ya: + /?' 



1- 

V'n-3 
Vn-3 

(px + qy) 

V n -4 



n-4 



(u n ) u (x,y) 
K-3 ) n . 

Un-l)22(y,X 



1 



/ rr(\ 



)(px + qy) 



u + v + (qx + py) 



V 



n-4 



n-4 



V r . 



n-l 



(60) 



The above equation reveals the recurrence relation among the partition functions of 
different Fibonacci generations from (n — 4) — th to (n + 2) — th.The partition functions 
have entered in the above equation through the quantities given by equations from (48) 
to (57) .This is in conformity with the symmetry properties of the Fibonacci chain. It is 
worth noticing that in the recurrence relation above ,the partition functions for both closed 
and open chains have entered. The procedure we have presented here through symmetry 
properties of the chain is generic to any substitution lattice. 

Illb Recurrence relation in the large N limit 

Recurrence relation (60) can be written in a much simpler form in the large iV limit. To 
achieve this we notice that in equation (60) both open and closed partition functions of a 
particular generation occur simultaneously.From the physical consideration that thermo- 
dynamic quantities must be same in the large iV limit for both open and closed partition 
functions of a particular generation ,we assume that the free energy per spin is same in 
both cases in a particular generation. Therefore, for the n — th generation relation between 
the partition functions can be written as: 



Z c n (F) = f n (r,f,H)Z° n (F) 



(61) 



where / n (r, r, H) is independent of the number of spins, N n + 1 ; N n being the number of 
bonds in the n — th generation. 

Case I: Ising Limit r = f ,H/0 

Considering equations (25) and (28) we can write equation (61) in the following form: 



MU-^ + ^(l-r)^ + 4(l-r)rTc^ifx A + 2 (A + 2 - A 2 ) ] Z ° 2m{F)U =" 

(62) 

which is equivalent to equation (15) in the large N limit where 2M stands for the number 
of spins in a large odd generation. From (61) and (62) we conclude that fn dd {r = f, H) is 
same for any large odd generation . Similarly, for large odd number of spins 2N + 1 in a 



10 



large even generation, equations (26) and (29) give: 



Z c 2N+l (F)\ e= - e =[Vt + 4(l-r)cosh(3Hx - 2 - ] (63) 



1 

A + (A + ~- A_ 



which is equivalent to equation (16) in the large N limit. From (61) and (63) we conclude 
that f^ ven (r = f, H) is same for any large even generation. 

Case II: r ^ r, H = 

Equations (25) and (28) give the exact expressions for open and closed partition func- 
tions with odd number of bonds in the following form: 

Z° 2M (F)\h=o = 2 2M {cosh (3ef L (cosh (3ef s (64) 

and 

Z C 2 M {F)\h=o = 2 2M [(cosh (3ef L+1 (cosh (3ef s + (sink f3ef L+1 (sinh f3ef s ] (65) 

It follows from equations (64) and (65) that in the large N limit 

Z C 2m(F)\h=o = (cosh Pe)Z° 2M (F)\ H=0 (66) 

Equations (26), (29) give the exact expressions for open and closed partition functions 
with even number of bonds as: 

Z° 2N+1 (F)\ H=0 = 2 2N+ \cosh Pef L (cosh (3e) Ns (67) 

and 

Z C 2N+1 {F)\ H=G = 2 2N+1 [(cosh (3ef L+1 (cosh (3ef s - (sink f3e) NL+ \sinh (3e) Ns \ (68) 

It follows from equations (67) and (68) that in the large N limit 

Z%n+i (F)\ H =o = {cosh (3e) Z2N+i{F)\h=o (69) 

From equations (61) , (66) and (69) we conclude that fn dd { r ^ f, H = 0) [fn Ven { r 7^ r, H = 
0)] is same for any large odd [even] generation . 

Considering the above cases we make the assumption that in the case of Fibonacci with 
H ^ 0, all fn Ven { r , t, H) in equation (61) for large odd generation number are equal 
and all f° dd (r,f,H) in equation (61) for large even generation number are equal; but 

f™ en (r,f, H)jtf° dd (r,f, H). 

Using the above properties it reveals that the expression Z°„ 1 Z^_ 3 — Z^_ ± Z°_ 3 appearing 
in equation (60) vanishes. As a consequence equation (60) becomes: 

v n-2\ , / . \ / \/„ / , ni v n- 



K+i)ii(?/, x)(a yx + Py X Ty—) + (u> n+1 ) 22 (y, x)(a' + P' yx Ty—) 

Vn-2 Vn-2 
1 V feven 

■ r^ n iw ~ , x (« + v + (qx + py)-^»)(^ - l)Z° n+l = (70) 
Vrr(l - i)(px + qy) V n - 2 V r 



where 
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Y y 1 D A 

(u n+1 ) n (y,x) = [ " " x x -=J-](y,x) x Z°_ 2 

1 n ly n ~ 1 n ly n Ty/T 



T n K' - T'K 



re ve il 
Jn+1 



n 



^(T n A' n - T> n A n/ 



{y,x)xz° n 



+1 



= &it- 2 (v,x)Z°_ 2 - $ 1 n 1 +1 (y,x)Z° n+1 (71) 
Here $(r, f, if)'s are indepedent of the partition functions. Similarly , one can write 

(u n+1 ) 22 (y,x) = ^_ 2 (y,x)Z° n _ 2 - <S>f +1 (y,x)Z° n+1 (72) 

In reducing (u n+ i) u (y, x) and (co n+ i) 22 (y, x) in the form (71) and (72) we have used the 
simplified expression 

fodd -I 
T Tl hi ± 

in _ r ( 7 o\ 

Vn - 1 

Since we have assumed (n — 2) — th generation to be even ;therefore (n — 1) — th , n — th 
generations will be odd and (n + 1) — th generation will be even. 
From equations (70), (71) and (72) we can write 

{^ +1 {y.x){a yx + /^fe) + ^ 2 +1 (y, x)(a' yx + /^fe) 

\rl seven 

+^(d)»[ U + v + (qx + py)^\{^t ~ l)}Z° n+l 
= {^_ 2 {y,x) + ^_ 2 {y,x)}Z°_ 2 



tf n+1 (r, f, )Z° +1 = M/ n _ 2 (r, f, tf)Z°_ 2 (74) 
In terms of free energy per spin equation (74) takes the form 

(N n+1 + l)F n+1 - (iV n _ 2 + l)F n+2 = (3 log^ (75) 

^ n+l 

Using the Fibonacci recurrence relation among the number of bonds expression (75) re- 
duces to 

(2r+l)F n+1 = F n _ 2 (76) 
where r = lim Nn ^^jj^ , the golden mean . 

Equation (76) shows that the free energy per spin for consecutive even generations are 
scaled by a factor (2r + 1) and therefore the thermodynamic quantities will be scaled in 
a similar fashion. 

IIIc Chemical Potential for H=0 

For calculating chemical potential we must know the exact form of the partition func- 
tion . Since we do not know the exact partition function for finite H we calculate it for 
H = . 

Let AT n _ 2 denote the number of bonds in the (n — 2) — th generation which is even . We 
add another spin at the end with a long bond , say , to the open chain. Then using 
equation (26) we get 

Nr Nq CT-i N L + 1 (J-, NS 

Z° Nn _ 2+2 (F)\ H=0 = r^f^Tr(l + ^) (1 + -j) (1 + a,) 

= 2 N ^ +2 (cosh pef L+1 (cosh (3ef s (77) 
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while the exact expression for the partition function with the open boundary condition 
in the (n — 2) — th generation is given in equation (67) as , 

Z° Nn _ 2+1 \ H=0 = 2 N ^ + \cosh (3e) NL (cosh (3ef s (78) 

For large generation the expression for chemical potential as defined in equation (17) is 

e^^)\ H=Q = v° ven (F)\ H=0 = fr 2+2 fl = teosh P<= 2 ^sh (3e) (79) 

The expression in the parenthesis stand for addition of a spin at the end of the chain with 
a short bond . 

In a similar fashion , the chemical potential for closed chain is 

v c even = Icosh /3e{or2cosh (3e) (80) 

We would have got the same result had we considered odd generation. Comparing expres- 
sions (79) and (80) with that of (22) we conclude that , in absence of external magnetic 
field the chemical potential do not depend upon the underlying lattice structure as long 
as we consider nearest neighbour interaction. 



IV Conclusion 

We have found an exact expression for the partition function of an Ising model on a 
reguler lattice with open boundary conditions in presence of magnetic field .This always 
includes closed partition functions because of the fact that out of four different spin 
configurations at the end points of the chain , two configurations |f and j j satisfy closed 
boundary conditions. Free energy per spin are equal for both open and closed chains in 
the large N limit , hence thermodynamic quantities are same in both cases ; but chemical 
potentials are different .For aperiodic chain we have established the recurrence relation 
among partition functions of different Fibonacci generations and we have also shown that 
this can only be obtained by using trace map relations and the characteristic symmetry 
property viz.," Mirror Symmetry" of Fibonacci generations . This procedure is generic to 
all substitution lattices. In our case the recurrence relation among the partition functions 
correctly reflects the fact that n — th and (n ± 6) — th generations are topologically same 
. As a consequence one must go through six times decimation renormalization group 
procedure to find scaling forms! of thermodynamic functions. The recurrence relation gets 
simplified in the large iV limit showing that free enrgy per spin in two consecutive large 
even generations are related through a scaling factor 2r + 1 . The expression for chemical 
potential in absence of magnetic field reveals that it is independent of the underlying 
lattice structure with nearest neighbour interaction . 
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